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Abstract 

We have studied hadronic r decay modes involving an rj meson using 490 fb^^ 
of data collected with the Belle detector at the KEKB asymmetric-energy e^e~ 
collider. The following branching fractions have been measured: B{t~ K~r]Ur) = 
(1.58 ± 0.05 ± 0.09) X lO^'', B{t~ K~-K^r]Ur) = (4.6 ± 1.1 ± 0.4) x 10~^ B{t~ 
TT~TT^r]Ur) = (1.35±0.03±0.07) X 10~^ ,B(r~ ^7r~K^r]Ur) = (4.4±0.7±0.2) x 10~^ 
and B{t^ K*~7]u^) = (1.34 ± 0.12 ± 0.09) x 10"^. These results are substantially 
more precise than previous measurements. The new measurements are compared 
with theoretical calculations based on the CVC hypothesis or the chiral perturbation 
theory. We also set upper limits on branching fractions for r decays into K~ K^rjUr, 
7r~ Kg-K'^rju-r, K^rjrjUr^ vr"?]?]!/,- and non-resonant K^T^^rjUr final states. 



1 Introduction 



Hadronic decays of the r lepton are very important for studying QCD phe- 
nomena at a low-energy scale. Various decay modes including an ?7-meson 
represent a wide class of such decays, which are still poorly studied. One of 
the important tasks in this field is to test effective chiral theories [1,2,3,4]. The 
r decay amplitudes in these theories are based on an effective Lagrangian of 
pseudoscalar, vector and axial- vector mesons with f/(3)L x f/(3)ij symmetry. 
On the other hand, for some decays the branching fractions as well as spec- 
tral functions of the final states can be predicted using the vector dominance 
model (VDM). Another important issue is testing the relations between the 
cross section of e^e" 7i^7i~r] and the spectral function in -K^TX^rjUj. 
decay predicted by conservation of vector current (CVC) and isospin sym- 
metry [5]. For all of the above tasks, precise measurements of the branching 
fractions as well as detailed studies of the final-state invariant mass distribu- 
tions are highly desirable. 
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Three r-decay modes involving r] mesons, r~ — K't^uJ^, t~ K*{892)~r]i'r, 
and T~ — s> tt'ti^tiVt-, have been studied by CLEO [6,7,8] and ALEPH [9]. 
However, the experimental statistics were very limited. For example, CLEO 
has the most precise measurements with 85 events in the r~ — > K'tjUt mode, 
125 events in the r~ Ti^-K^riVr mode and 25 in the r~ — * K*~'r]Ur mode. 
The limited statistics of these measurements result in large uncertainties for 
the branching fractions and do not allow study of the hadronic mass spectra. 

We report a detailed study of the branching fractions for K'tjUt-, 
K'TT^TjUr, Ti'TT^rjVri as Well as Ti'K^rjVr decays. The K'tt^ invariant mass 
in the K~Ti^riVT- decay and the i^'K^ invariant mass in the n'KgTjUr decay 
are studied to evaluate the branching fraction of the r~ — >■ K*{892)~'r]h'r decay. 
In addition, some other final-state invariant mass distributions are presented 
and discussed. 

This study was performed at the KEKB asymmetric-energy e~^e~ collider [10] 
with the Belle detector. The results are based on a 490 fb~^ data sample that 
contains about 450 million r-pairs, which is almost 100 times larger than in 
any of the previous measurements. 

The structure of this paper is as follows. Section 2 describes the Belle de- 
tector. Event selection is explained in Section 3. The branching fractions are 
determined in Section 4, where that of r~ A'* (892) "771/,- is evaluated from 
analysis of the K^-K^rjUT- and Ti'K^grjVj. samples in Section 4.3. Finally, 

some concluding remarks are given in Section 5. 



2 Experimental apparatus and Monte-Carlo simulation 

The Belle detector is a large-solid-angle magnetic spectrometer that consists 
of a silicon vertex detector (SVD), a 50-layer central drift chamber (CDC), an 
array of aerogel threshold Cherenkov counters (ACC), a barrel-like arrange- 
ment of time-of- flight scintillation counters (TOF), and an electromagnetic 
calorimeter (ECL) comprised of CsI(Tl) crystals located inside a supercon- 
ducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- return 
located outside the coil is instrumented to detect mesons and to identify 
muons (KLM). The detector is described in detail elsewhere [11]. Two inner 
detector configurations were used. A 2.0 cm radius beam pipe and a 3-layer 
silicon vertex detector were used for the first sample of 140 fb~^, while a 1.5 
cm radius beam pipe, a 4-layer silicon detector and a small-cell inner drift 
chamber were used to record the remaining 350 fb~^ [12]. 



Unless specified otherwise, charge-conjugate decays are imphed throughout the 
paper. 
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Particle identification (PID) of the cliarged tracks uses likeliliood ratios, Vx-, 
for a charged particle of species x (x = /i, e, K, n or p). Vx is defined as 
Vx = Lx/T,yLy (the sum runs over the relevant particle species), where Lx 
is a likelihood based on the energy deposit and shower shape in the ECL, 
the momentum and dE/dX measured in the CDC, the particle range in the 
KLM, the light yield in the ACC, and the particle's time-of-fiight from the 
TOP counter [13]. The efficiencies for kaon, muon and electron identification 
are 77% for momenta of p > 0.3 GeV/c, 92% for p > 1.0 GeV/c and 94% for 
p > 0.5 GeV/c, respectively. 

The detection efficiency for each signal mode and the level of background (BG) 
contribution are estimated from Monte Carlo (MC) simulations. To generate 
signal events as well as r-decay backgrounds, the KKMC program [14] is used. 
Por the decay modes not covered by KKMC, like r~ — > K'TjUr, K*{892)~rii'r, 
. . ., the signal MC is produced assuming a pure phase-space distribution of the 
hadronic system with a.V — A weak interaction. The background from e~^e~ 
qq {qq continuum) is simulated using the EvtGen procedure [15]. The detector 
response is simulated by a GEANT3 [16] based program. Throughout the 
paper the efficiencies quoted include the branching ratios of the corresponding 
r] decays. 



3 Event Selections 

The signal events should have the following common features in the e'^e" 
''"signal ''"tag reactiou: 

^signal ^ X~ + 7] + n{< 1)7 + (missing), 
nig ^ We)^ + n(< 1)7 + (missing), 

where X' denotes K' , K'n^, 7r"7r°, 7r~K'^g or K*(892)- systems. The 7r°, 
Kg and K*{892)~ are reconstructed through tt*^ — > 77, Kg — > tt^tt", and 
K*{892)~ Kg7r~ , K'n^ decays, respectively. The 77 meson is identified 
through the 77 — > 77 decay mode for all r channels while the 77 — > 7r~^7r~7r° 
mode is added for the r~ K'rju decay mode. 

The signal events, therefore, comprise either two or four charged tracks with 
zero net charge and two or four 7's. In order to take into account initial-state 
radiation, we allow at most one extra 7 on both the signal and tag sides. Each 
charged track is required to have a transverse momentum pt > 0.1 GeV/c, and 
a polar angle of —0.866 < cos 6 < 0.956, where pt and 6 are measured relative 
to the direction opposite to that of the incident beam in the laboratory 
frame. Each photon candidate should have an energy i?^ > 0.05 GeV within 
the same polar angle region as for the charged particles. 
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The tag-side rtag is required to decay into leptons, i.e. r"*" E'^vgi'r = 
e/fi), which corresponds to a branching fraction of 35.2% [17]. The lepton 
candidates are required to be well identified by requirements on the Ve or "P^ 
PID parameters with p > 0.7 GeV/c. 

The thrust axis (rithrust) is defined to maximize Khmst = Si|pf^-nthrust|/Si|pf^ 
where pf ^ denotes the momentum of each particle in the center-of-mass sys- 
tem (CM), and the sum runs over all detected particles in an event. The 
requirement Khrust > 0.8 is applied to remove qq events. An event is divided 
into two hemispheres by the plane perpendicular to rithrust- The hemisphere 
that includes Tgignai with an r] is referred to as the signal side, while the opposite 
hemisphere, which includes rtag, is defined as the tag side. 

In order to remove non-r pair backgrounds, such as the qq continuum, the 
total energy in CM is required to be in the range 3.0 GeV < -Et^tai < 10.0 
GeV. In addition, the invariant masses of both the signal and the tag sides are 
required to be smaller than the r-mass: Msig, Mtag < rrir (1.78 GeV/c^). To 
remove two-photon backgrounds, the missing momentum should correspond to 
a particle crossing the fiducial region; we require that the missing momentum 
satisfy —0.866 < cos6'(pmiss) < 0.956. After applying these requirements, the 
qq background level is 3% of the signal yield. Background from two-photon 
reactions is negligible. 

To reconstruct an rj, two 7's with > 0.2 GeV are required in the barrel 
region (—0.624 < cos 6* < 0.833). In addition, we allow combinations with 
at most one extra 7 with 0.05 < < 0.2 GeV in the endcap region on 
the signal side to take into account initial-state radiation as well as beam- 
induced BG clusters in the calorimeter. In order to reduce the number of 
incorrect combinations with a 7 from decay (denoted hereafter as 7,^0), the 
r^-candidate 7 (7^) should not form a 7r° mass with any other 7, i.e., a vr^-veto 
is applied. The tt^ mass window is defined in this paper as 0.105 GeV/c^ < 
< 0.165 GeV/c^, which is a ±3(T range of the detector resolution. 

3.1 T~ —>■ K'TjUr decay 

Candidate events for this decay mode must contain one charged track and at 
least two 7's in the 77 — > 77 case, or three charged tracks and two 7's in the 
r] — >■ -7r"''7r~-7r° (vr*^ — > 77) case. 

A kaon is identified by the Vk parameter with p > 0.3 GeV/c. In addition, 
a low Ve value is required to remove beam electrons from two-photon pro- 
cesses; this requirement reduces two-photon BG to a negligible level. The 
opening angle between the K~ and r] is required to satisfy the condition 
cos6'(P^'^; Pi"^) > 0.8, to reduce the combinatorial BG. In addition, the 
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Fig. 1. (a) and (b) M^+^-^o distributions for K'rjv-r- selection in 77 —> 77 and 
T] — > 'k'~'k~^it^ decays, respectively. Data are fit with a Crystal Ball function plus 
a second-order polynomial for the BG. The best-fit result is indicated by the solid 
curve with the BG shown by the dashed curve. 



opening angle and energy of the two 7^'s should fulfill the following condi- 
tions: 0.5 < COS^(pCM.pCM) < Q_9g I^CM _ ^CM|/|^CM ^ ^CM| < g.S, 

respectively. 



For rj — s> n^Ti^TX^ reconstruction, the selection criteria are basically the same 
as in the 77 — 77 case. Here 77 candidates should have two additional charged 
tracks, and two of the 7's must form a vr'' instead of an rj. 



The 77 and ti^ti^tx^ invariant mass distributions around rj mass are shown 
in Figs.l (a) and (b), respectively. The rj peak is clearly seen in both cases. 
The selection efficiencies (including the intermediate branching fractions of 
^{f] 77), -6(77 7r+7r"7r°), and B{t~ I'v^v^)) are e = 0.94% in the 
?7 — i> 77 and e = 0.16% in the 77 tt+tt^tt^ case. 
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likelihood ratio likelihood ratio 

Fig. 2. Distribution of the likehhood ratio for (a) K~TT^r}UT and (b) n'TT^'qi'T-. Dots 
are experimental data. The hatched and normal histograms indicate the signal and 
BG MC distributions, respectively. The dashed line shows the threshold of this 
selection. 

3.2 K'n^Tju^ and TX'Tx^rju^ decays 



For these decays, an 77 is reconstructed through the 77 mode only, so that the 
signal event should contain one charged track and four 7's on the signal side. 

All 7's, except for a possible extra 7, are required to be detected in the 
barrel region of the calorimeter, and two of them should form a combina- 
tion consistent the 7r° mass. The total momentum on the signal side is re- 
quired to satisfy Z^pSg' > 2.5 GeV/c. A condition on the opening angle be- 
tween the missing momentum in the CM and the thrust axis is imposed: 
cos6'(p^^g; n^jfj.ygj < —0.6, to remove BG's with large missing momentum on 
the signal side. Charged kaons and pions are selected by requiring large Vk 
and small Ve in the K'Ti^rjVj. mode, while small Vk and small Ve are selected 
in the tt'tt^itiVj. mode. 

For further background rejection, we construct a likelihood using the following 
seven variables: Vthrust, momenta of the 77 and vr" in the CM (pj^^, p^o^), the 
missing-mass squared (M^j^g), the energy of the 7^'s in the CM [E^^), Y.P^^^ 
and coiiOijfj^q^^p^^). The MC simulation is used to study the likelihood dis- 
tributions of generic r+r~ and the small residual qq BG {L-qq) as well as the 
signal (i^sig)- A likelihood ratio is defined as R = Lsig/(Lsig + -^bg)- The R 
distribution is shown in Figs. 2 (a) for K'tt^tiVj. and (b) for tt^vt^?]!/,- decay; 
with the requirement R > 0.6, about half of the background is removed, while 
93% and 90% of the signal samples are retained, respectively. 

The resulting M^^ distributions are shown in Figs. 3 (a) and (b) for the K'TT^rjVj. 
and 7r~7r°?7i/^ modes, respectively. The 77 peaks are clearly visible. The detec- 
tion efficiencies are e = 0.35% and 0.47% in the r~ — > i^~7r°?7z/^ and n~7i^r]h'r 
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Fig. 3. M^^ distributions for (a) K~TT^r]UT- and (b) -k^t^^i^Vt candidates. Data are fit 
with a Crystal Ball function plus a second-order polynomial for the BG. The result 
of the best fit is indicated by the solid curve with the BG shown by the dashed 



curve. 



decays, respectively. 



3.3 -K'K^siqVr decay 



On the signal side, three charged tracks and at least two 7's are required. 
To suppress two-photon BG, the tracks should not be identified as electrons. 
Without K/ti id the pion mass of m^± = 0.13957 GeV/c^ is assigned to all 
three tracks on the signal side, two of the three tracks with opposite charges 
are the daughters, the invariant mass (M^+^-) of which is required to be 
within a mass window (3(t), 0.4907 GeV/c^ < M^+^- < 0.5047 GeV/c^. To 
ensure that these tracks are daughters, the distance from the interaction 
point to the decay vertex on the plane perpendicular to the beam axis 
('"vertex) is required to satisfy the condition: 0.5 cm < rvertex < 30 cm. To 
avoid multiple candidates, the mass of the other possible combination must 
not satisfy the mass condition. Pairs of 7's with > 0.3 GeV are used 
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to form 7] candidates. We require that the energy of an r] candidate, E^j, must 
be greater than 1.2 GeV. 

The resulting M^^ distribution is shown in Fig. 5 (a). The detection efficiency 
is e = 0.32%, which includes the tag-side efficiency, B{rj — > 77) and B{Kg 



4 Branching Fractions 



4-1 K rjVr, K -K^rjVr and TT Tx^rjUr decays 



4.1.1 Signal candidates 

The rj meson yields are obtained from fits to the 77 invariant mass distri- 
butions shown in Figs.l and 3. Fits are performed with a Crystal Ball (CB) 
function [18] plus a second-order polynomial function to represent the 77 con- 
tribution and combinatorial BG, respectively. All five parameters of the CB 
function, as well as three others for BG are treated as free parameters in the 
fit. 

The results of the fit are shown by the solid curve in the corresponding fig- 
ures. The best fits give r] masses of = 0.5449 ± 0.0006, 0.5444 ± 0.0010 
and 0.5443 ± 0.0004 GeV/c^ for K'tjUr, K'tt^tjUt, and n'lr^riUr decays, re- 
spectively, in the 77 — 77 subsample, while = 0.5474 ± 0.0007 GeV/c^ for 
K'TjVr in the 77 — > 7r+7r~7r° subsample. The corresponding mass resolutions are 
a^^ = 0.0110 ± 0.0006, 0.0107 ± 0.0015, 0.0116 ± 0.0003, and 0.0075 ± 0.0004 
GeV/c^, respectively. These values of and cta//^ are in good agreement with 
those obtained from the MC simulation, although the obtained masses are 
shifted by 2 - 3 0" from the 77 mass in Ref. [17], which is due to the incomplete 
detector calibration. For instance, M^^ = 0.5442, 0.5450, 0.5464, and 0.5473 
GeV/c^ for a generated mass of = 0.54745 GeV/c^. 

The rj yields obtained from the fits are N^j = 1,387 ± 43, 270 ± 33, and 
5, 959±105 events for the K'rju^, K'Tr^rjUT-, and -K'Tr^rjVj. modes, respectively, 
in the i] ^ ^'-f case, while A'^^ = 241 ± 21 events for the K'-rjUr decay in the 
77 — >• Tr+vr"??^ case. These yields include BG, attributed to other r decays with 
?7 meson(s), such as the r~ — > 7r~7r°7r°7/i/^ and n'lr^n'riVr] QQ continuum with 
T) meson(s); cross-feeds among the signal modes. 
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4- 1-2 Background 



The possible backgrounds in the K'-rjUr, K~7r rjUr and tx^tx rjVr decay modes 
include (a) feed-across from signal modes and (b) other processes such as 
T — i> TT^TT^ir^rii'r, K'rjrjVr-, T^'W^T, Ti^riVr and qq. We first discuss the second 
source. 

BG from generic r'^r" decays is evaluated using a MC simulation where the 
value of B{t^ tt^tt^tt'^?]//,-) is taken from Ref. [19]. This BG is negligible for 
K'tjUt- and K^-r^itiVt- decays, and is only 72 ±20 events, which is smaller than 
the statistical uncertainty in the rj yield for the ti'ti^tiVt decay. 

To examine the r~ K^rjrjVT- and 'K"ririVj. contaminations, we analyze data 
using selection criteria similar to those used for the K^ir^rjUT- and ir^-K^'qUj. 
modes. An additional r], instead of a 7r°, is reconstructed, and no selection on 
the likelihood ratio is implemented. Requiring one photon pair to be in the rj 
mass region, 0.48GeV/c2 < M^l} < 0.58 GeV/c^, which includes the largest 
energy 7, the invariant mass distribution of the other 77 candidate, (M^^^), is 
shown in Figure 4. The evaluated 77 yields are lAtli and 4.4 ± 6.8 events 
in the K'Tj-rjUr and ir'TirjUr decays, respectively, with detection efficiencies of 
e = 0.19% and 0.22%, respectively. We then set the upper limits [20] 

B{t~ R-rjrjUr) < 3.0 x 10~^ (1) 

and 

B{t^ Ti~r]r]Ur) < 7.4 x lO"'' (2) 

at the 90% confidence level (CL) including systematic uncertainties of 6.9% 
and 7.1%, respectively. Since the cross feed probability for these modes with 
extra ?7's is below 3%, their contamination is ignored. 

The r~ tt^tjUj. decay proceeds via a second-class current so that its branch- 
ing fraction is expected to be small, as low as 10~^ [3] or less. Since the 
efficiency for detecting n'TjUr as signal is of the order of 10~^, the total con- 
tamination is then ~ 10~^ or less so that this decay can be ignored. 

BG from the qq continuum is examined using both MC and data in order to 
take into account the 77 production uncertainty of the qq MC. A qq enriched 
sample is produced with some variations of the signal selection criteria. We 
require M^^g > Tn^-, while the requirement Mgig < rrij. is not applied; the PID 
requirement on the tag side is reversed. Other criteria are not changed. A 
comparison between the resulting rj yields of MC and the data resulted in 
the following MC scale factors: 2.3 ± 0.7 for K'rjVri 2.6 ± 0.5 for K'Tr'^rjUr, 
and 2.9 ±0.3 for 7r"-7r°?7z/^. Therefore, we obtain rescaled qq contaminations of 
39.7 ± 15.8 events in K'Tju^-, 212 ± 29 events in tt'tt^i^Pt and 27.0 ± 8.5 events 
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Fig. 4. Mass of the 7,j's combination, which does not include the largest energy 7,,, 
M-^-y for (a) the t~ K~rir]Vr sample and (b) the t~ 7r~r]r]i'-j- sample. Data 
are fit with a Crystal Ball function plus a second-order polynomial for the BG. The 
result of the best fit is indicated by the solid curve with the BG shown by the dashed 
curve. 

in the K'n^rjUr decay. These BG's correspond to 2.8%, 10% and 3.6% of the 
raw signal yields, A^^, respectively. 

The yields in the three candidate samples after the BG subtraction evaluated 
above, Ni, {i = K~iriv, K^ir^rju, TT^ir^rju) are used to evaluate the correspond- 
ing branching fractions, Bi. To take into account the cross- feeds between decay 
channels, we solve the following system of linear equations: 

iV, = 2iV..^e;.i3, (j = 1,2,3), (3) 

i=l 

where e*- is the efficiency for detecting mode j as mode i, which is estimated 
by MC; Nrr is the total number of r-pairs produced. In the 77 — 7r"'"7r~7r'^ case, 
B{t~ K^rjUr) is calculated using the relation N[ = 2Ar,, Eti e'l^S, (the 
prime indicates quantities for this 77 mode) with B2 and B3 obtained in the 
?7 — *■ 77 case. The statistical uncertainty in A''( is ~ 10%, while the cross-feed 
contamination is ~ 4%, and is thus insignificant. 

The resulting branching fractions are B{t~ K'rjVr) = (1.57 ±0.05) x lO"'^, 
B{t- K-Ti^r]P^) = (4.6 ± 1.1) X 10"^ and B{t^ Tr-yrVr) = (1-35 ± 
0.03) X 10-3 in the ^ 77 case, and B{t^ K^Wr) = (1-60 ± 0.15) x 10"^ 
in the rj — ^ 7r"'"7r-7r'' case. 

4- 1-3 Systematic uncertainties 

Systematic uncertainties are discussed below, and listed in Table 1. 
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Table 1 

Summary of systematic uncertainties in each mode (%). The uncertainty in the 
number of t~^t~ events comes from the uncertainties in our luminosity measurement 
by Bhabha events and the cross section of e^e~ t'^t^ [21]. 



Signal modes 


K r]u 


K 'K^l]U 


TT ■K^1]U 


K r]u 


Items 






BG subtraction 















0.6 


1.8 X 10-3 


- 




0.3 


_ 


4.2 X 10-2 


0.4 




7.5 X 10~2 


3.3 


_ 


0.1 










0.4 


- 


qq 


1.5 


6.0 


0.5 


1.5 


Detection efhciency 










K/iT- /lepton-id 


3.3/ 2.3 


2.2/ 2.8 


1.0/2.6 


2.8/ 2.6 


Tracking 


1.3 


1.3 


1.3 


3.3 


TT^ /r] —>■ 77 


-/ 2.0 


2.0/ 2.0 


2.0/ 2.0 


2.0/ - 


TT^ veto 


2.8 


2.8 


2.8 




Stat, error of signal MC 


0.5 


1.7 


0.5 


1.3 


Bit] vr+^-^O) 








1.6 


Luminosity meas. 

a{e~^e~ — > r+r^) 


1.4 
0.3 


Total 


5.9 


9.1 


5.3 


6.2 



Subtraction of the BG discussed above provides the systematic uncertainties 
includes the statistical errors in the detection efficiencies for BG processes, as 
listed in the Table. Sizable uncertainties are found only in K~iv^rjv decay from 
Ti'Ti^rjVr and the qq continuum and amount to 3.3% and 6.0%, respectively. 

The systematic uncertainties for PID, track finding, 7r° and rj reconstruction, 
and the 7r°-veto are also estimated. The modeling of the PID likelihood is 
tested by studying inclusive D*~ samples for K/tt and two-photon 77 — > t^l~ 
samples for leptons; the Uncertainty in the PID efficiency is 2-3%, 1% and 
~ 2.5% for K^, TT^ and {fi/e)"^, respectively. The uncertainty for finding a 
track is 1.0 % for hadrons and 0.3% for leptons. The uncertainty in ti^/t] '-ff 
reconstruction is 2.0%, obtained from a comparison of 77 — 77 and rj vr^vr^vr^ 
data samples. The vr^'-veto inefficiency in tt'-k'^tjUt- samples is evaluated by 
comparing its effect on the data and MC samples when the tag-side lepton- 
PID criterion is reversed from the ordinary selection. 
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The statistical errors on the MC simulation are 0.5%, 1.7% and 0.5% for 
T~ — s> K'tjUt-, K'TT^riVj. and n'-K^rjUj. in the rj ^ 'y'y case, respectively, and 
1.3% for T~ — > K'TjUj. in the rj — s> 7r+7r~7r° case. The branching fractions 
used in the MC are cited from Ref. [17], where B{ri — > Tr'^Tr^vr'') yields a 
sizable uncertainty of 1.6%, while others are negligible. The uncertainty of the 
a{e^e~ t^t~) cross section and the integrated luminosity are 0.3% and 
1.4%, respectively. 

All contributions are summed up in quadrature to obtain the total uncer- 
tainties; they amount to 5.9%, 9.1% and 5.3% for B{t~ K'tiUt-), B{t~ 
K'TT^rjVT-) and B{t~ rr'n^riu^), respectively, in the rj ^ '-ff case, and 6.2% 
for B{t~ K'tjUt-) in the rj tt+tt'tt'^ case. 

Taking into account the systematic errors, we obtain the following branching 
fractions: 

= (1.57 ±0.05 ±0.09) X 10"^ for r/ ^ 77, (4) 
= (1.60 ±0.15 ±0.10) X 10"^ for r/ ^ 7r+7r"7r°, (5) 

= (4.6 ± 1.1 ±0.4) X 10-^ (6) 

= (1.35 ± 0.03 ± 0.07) X 10"^ (7) 
For K'rjUr, the two measurements are combined to obtain 

= (1.58 ±0.05 ±0.09) X 10"^ (8) 



4.2 IT KgrjUT- decay 



4-2.1 Signal candidates 

For 7r~Kgr]i'T-, the signal yield A^^ is evaluated from a fit to the M^^ mass 
distribution in the same way, as discussed in 4.1.1, but the parameters of the 
CB function are fixed to those determined by the signal MC simulation. The 
fit result is shown in Fig. 5 (a), yielding = 161 ± 18 events. 

We examine r~ — K'KgfjUr and tt^tt+tt"?]!/^ decays in order to determine 
their cross-feed contaminations, with selection criteria similar to those for 
Ti~ KgTiVr- For K~ KgrjUT, the charged track on the signal side is required to 
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0.5 0.6 0.7 0.4 0.5 0.6 0.7 



Myy(GeV/c ) 



Myy(GeV/c^) 



Fig. 5. M^^ distributions for (a) tt~ KgrjVT-, (b) K~ K^rjUr, (c) tt~ tt^ tt~ tj^t without 
a Kg, and (d) tt' tt^ Kgrju-r samples. Data are fit with a Crystal Ball function plus 
a second-order polynomial for the BG. The result of the best fit is indicated by the 
solid curve with the BG shown by the dashed curve. 



be a kaon; A'^ = 5.6 ± 3.6 events is obtained from the fit (see Fig. 5 (b)). 
For T~ — > 7r~ 7r~^ 7r~ rjUT- , its M^^ mass distribution, as shown in Fig. 5 (c), is 
formed from a sample in which the tt+tt" invariant mass lies outside the Kg 
mass range while the other requirements for Kg daughters are the same; the 
fit gives iV^ = 67.9 ± 17.0 events. 



4-2.2 Background 

The possible background to the signal comes from 7r~ ir^ Kg-qUr and the qq 
continuum. The former sample is selected, requiring an additional vr'' on the 
signal side. No clear rj peak is found in the M^y distribution, as can be seen 
in Fig. 5 (d): N^j = 4.715;? events is obtained with a detection efficiency of 
e = 0.08%. We therefore set an upper limit on the branching fraction B{t~ — > 
Tr-K'^ir^rjUr) < 2.5 x 10"^ at the 90% CL. The BG from this mode in the 
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TT Kgrju-j- sample is calculated to be 1.6Ii;7 events. 

The contamination of qq in ti' K^rju-r decay is 21.1 ±7.5 events, following from 
MC simulation with a scale factor estimated from the qq enriched sample. 

After subtracting the two above BG's from the yield A^,, for ti' K^rjVr decay, 
we solve the simultaneous equations of Eq.(3) in order to take into account the 
cross-feeds among the three modes. Now, the (super)subscripts i, j = 1, 2, 3 for 
the equations correspond to 1 = n'Kg'qi'T-, 2 = K~ Kgrju^-, 3 = Tr^Tr+vr^r/z/T-. 
The last category includes other peaking BG from generic r decays, which are 
estimated using samples in the Kg sidebands. The detection efficiencies are, 
for instance, e\ = 0.32% and = 0.23%. The resulting branching fractions 
are B{t- n-K^-qPr) = (4.38 ± 0.75) x 10"^ and B{t~ K-K^riVr) = 
(1.9 ±2.0) X 10"^ respectively. 

4-2.3 Systematic uncertainties 

Systematic errors for B{t~ TT'K^rjVr) include the following sources. Un- 
certainty in the signal CB function used to fit the M^-^ spectrum is examined 
by varying and ctm^ within their errors: the fitted yield A^^ varies by 0.7%. 
The dominant BG comes from the qq continuum: its uncertainty is evaluated 
to be 6.2%, while that of the others is negligible. The uncertainties in the 
luminosity evaluation and cr(e"'"e~ t^t~) are 0.3% and 1.4%, respectively. 
The value of B{K^ — > tt^tt') used in MC has a 0.4% error [17]. The uncer- 
tainty in the CB function yields a 0.5% contribution. The total systematic 
uncertainty is consequently calculated to be 7.9%, by adding all of the above 
errors in quadrature. 

As a result, we obtain the following branching fractions: 

B{t~ ^ Ti^Klr^Ur) 

= (4.4±0.7±0.3) X 10-^ (9) 

and 

B{t~ K-Kl7]Vr) < 4.5 X 10-^ (10) 

at the 90% CL. 

4-3 T~ K*~r]VT- decay 



The [Kit) invariant mass spectra in both r — A' ir^rjUr and r tt KgrjUr 
decays are analyzed to determine the A'*(892)~ content in the final state. 
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0.8 1.0 

Mk„o (GeV/c^) 




Mk,„ (GeV/c2) 

Fig. 6. (a) K~TT^ invariant mass distribution for K~iT^rjVr events and (b) the vr~i^^ 
invariant mass distribution for tt~ K^i]Ur events. Data are fitted with a convoluted 
Breit-Wigner function, and the best fit is indicated by the sohd curve, while the BG 
component is shown by the dashed curve. The fits give Nx*-nu.r = 122 it 17 events 
(a) and Nj^,-^^^ = 123 ± 14 (b), respectively. 

4^.3.1 — > K^Tx^TjUT samples 

To select r~ K'TT^rjv^, the rj signal is identified with the criterion 0.50 GeV/c^ < 
My^ < 0.58 GeV/c^, and then the K~7r^ invariant mass M^-nO distribution is 
analyzed (see Fig. 6 (a)). The contribution of the combinatorial BG under the rj 
peak is evaluated from the sidebands with 0.43 GeV/c^ < M^^ < 0.48 GeV/c^ 
and 0.60 GeV/c^ < M^^ < 0.65 GeV/c^. The dashed curve in Fig. 6 (a) indi- 
cates the BG component estimated from the above samples and from MC that 
includes a r decay with an rj. A MC study indicates that the enhancement 
in the higher M^-^o mass region mostly arises from tt'tt^itiVj. events through 
misidentification of vr" — s> K~ , while the dominant BG, mostly populating the 
lower mass region, is composed of contributions from generic r decays. A clear 
-ft'*(892)~ signal over BG can be seen in the spectrum shown in Fig. 6 (a). 

The observed Mx-ttO spectrum is fitted with a K* Breit-Wigner (BW) function 
plus a BG component, where the former is convoluted with a response function 
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with a mass resolution of o"m^_^o — 13.5 MeV/c^. The latter is fixed so as to 
reproduce the M^-^o spectrum of the sidebands. The mass and width of the 
BW function are fixed to those of the K*{892)- [17] in the fit. The best-fit 
result is indicated by the solid curve in the figure, and gives Nx*-nur = 122±17 
events with a xVd-o.f = 27.4/27. 

We examine the non-resonant K~n^ contamination by adding an additional 
term in the fit, assuming a phase-space distribution with a.V — A weak in- 
teraction for a hadronic final system. The fit gives Nx—r^v-r = 104 ± 20 events 
and A^non-K*- = 42±25 events with a x^/d.o.f = 24.5/26. Since no significant 
difference between the two fits is found within the errors, we here simply take 
the intermediate state to be purely K*~riVj. and the difference is taken into 
account as the systematic uncertainty for the fitting function. Even allowing 
for interference between the K*~ and non-resonant /r~7r°, the difference is 
found to be negligible. 

No significant i^*(892)~ BG contribution is found in generic r decays, while 
6.5 ± 2.3 events are estimated from the qq MC. 

The detection efficiency is evaluated from MC as e = 0.12%, including the 
branching fractions for i^*(892)~ — > K~ti^ and t] ^ 'j'y decays. 

The systematic uncertainties in the evaluation of B{K* {892)~T]h'T-) are summa- 
rized in Table 2. The dominant error arises from the uncertainty in the fitting 
function discussed above. It amounts to 15%. Other sources of systematic un- 
certainties are the same as those for K~7r'^rji>T- decay. The total systematic 
uncertainty is 16.2%. 

Thus we obtain the branching fraction 

= (1.13 ± 0.17 ± 0.18) X 10"^ (11) 

For the non-resonant K'n^rjVT- decay, we set an upper limit on its branching 
fraction 

B{t~ A'-7r°r/Z/,)non-rcsonant < 3.5 X 10"^ (12) 

at the 90% CL, assuming a pure phase-space distribution for the final hadronic 
system. 

4-3.2 T~ ^ K^g'K~riVj. samples 

The Mj^o^^- distribution is shown in Fig. 6(b) and is similar to the M^-^o 
case. A clear ii'*(892)~ signal can be seen over a small continuum BG. 
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Table 2 

Summary of the systematic uncertainties in K*'~rjVr analysis (%). 



A (892) decay mode 


T^— 

K vr^ 




BG subtraction 
qq 


2.0 





Detection efficiency 






K/tt- /lepton-id 


2.2/ 2.5 


-/ 2.4 


Tracking 


1.3 


3.3 




2.0/ 2.0 


-/ 2.0 


TT^'-veto 


2.8 




Stat, error of signal MC 


1.7 


0.7 






0.1 


Mass spectrum 


0.5 


0.5 


Luminosity meas. 


1.4 






0.3 




Fitting function 


15.0 


4.1 


Total 


16.2 


6.4 



Assuming no non-resonant KgTT~ contribution, the distribution is fitted with a 
convoluted BW function plus a BG, composed of a third-order polynomial, ap- 
plying the same method used for the K^-k^tiVt- sample. The yield of i^*(892)~ 
is 123 ± 14 events. The K*~ BG is estimated to be 14.5 ± 5.0 events from the 
qq continuum, and is subtracted from the above yield. 

We examine a possible non-resonant A'^vr" contribution, whose mass spectrum 
is calculated by MC, assuming a pure hadronic phase-space distribution with 
a.V — A weak interaction. Including the non-resonant component, the best fit 
gives Nk*- = 121 ± 16 events and A^non-resonant = 3 ± 15 events. Therefore, we 
give results assuming no non-resonant background. 

The detection efficiency is evaluated by MC as e = 0.10%, which includes the 
relevant branching fractions: B{K* (892)" n^Kg), B{Kg tt+tt"), and 
B(r] 77). 

The systematic uncertainties are summarized in Table 2. Their magnitudes are 
similar to those in the K'n'^rjUr case, except for the non- resonant contribution. 
The total systematic uncertainty is 6.4%. Consequently, the branching fraction 



19 



Table 3 

Comparison with previous results. 





Branching fraction ;B (xlO ^) 


Mode 


This work 


Previous exp. 


Reference 




1.58 ±0.05 ±0.09 


2.6 ±0.5 ±0.5 


CLEO [7] 






2.9 ± 1.3 ±0.7 


ALEPH [9] 


T IT -K rjVr 


13.5 ±0.3 ±0.7 


17±2±2 


CLEO [6] 






18±4±2 


ALEPH [9] 


T~ — > K~'K^r]i'j- 


0.46 ±0.11 ±0.04 


1.77 ±0.56 ±0.71 


CLEO [8] 




0.44 ±0.07 ±0.02 


1.10 ±0.35 ±0.11 


CLEO [8] 


T~ K*~r]l'r 


1.34 ±0.12 ±0.09 


2.90 ±0.80 ±0.42 


CLEO [8] 



is 

= (1.46 ± 0.16 ± 0.09) X 10"^ (13) 

Two measurements using K*{892)~ —> K~7r^ and KgTT~ decays are in agree- 
ment, in accordance with isospin symmetry. Therefore we combine the two 
results and obtain 

B{t- K*{m2)-7]Vr) 

= (1.34 ± 0.12 ± 0.09) X 10"^ (14) 



5 Result 

Using a high statistics 450 million r-pair data sample from Belle, we have 
obtained the following branching fractions for five different decay modes: 

B{t- K-r]u^) = (1.58 ± 0.05 ± 0.09) x 10-^ 

B{t^ TT^TrVr) = (1-35 ± 0.03 ± 0.07) x 10-^ 

B{t~ A'^TrVr) = (4.6 ± 1.1 ± 0.4) x 10~^ 

B{t~ tx~KIt]u^) = (4.4 ± 0.7 ± 0.3) x 10-^ 

B{t- K*{m2)-r]Ur) = (1.34 ± 0.12 ± 0.09) x 10-^ 

where the first and second errors are statistical and systematic, respectively. 
We also set the upper limits on the following decay modes at the 90% CL: 

B{t~ K^K^r^Vr) < 4.5 x 10~^ 
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B{t- Ti-KlTi^r]Vr) < 2.5 x 10'^, 

B{t- K~T]r]Ur) < 3.0 X 10"^ 

B{t- 7r-rir]Ur) < 7.4 x 10'^ 

B{t~ /'i^~Vr°?7Z/^)non-resonant < 3.5 X 10 



In Table 3, our results are compared to those previously obtained by the 
CLEO [6,7,8] and ALEPH [9] collaborations. It is clearly seen that the preci- 
sion of the measured values has been considerably improved. 

It is also noteworthy that the central values of our branching fractions are 
in all modes lower than those of the other experiments [6,7,8,9]. This fact 
can be mostly attributed to the underestimation of the BG contamination in 
previous low statistics measurements. For instance, the qq BG estimation in 
the previous analyses relied on the MC. However, this analysis evaluates the 
BG rate using a qq enriched data sample. Furthermore, Ref. [7] ignored the 
K~7r^r]Ur BG contribution, since they estimated it as ~ 0.1% BG, based on a 
theoretical calculation that predicted its branching fraction to be 8.8 x 10~^ 
[3]. However, its rate is about 1% of the signal in our measurement, and its 
branching fraction is about five times as large. Our high statistics allows us to 
reliably and precisely scrutinize various BG contributions, while the previous 
measurements, which had low statistics, did not have sufficient sensitivity to 
notice, or correctly judge, the BG contamination. 

Our results are compared with different theoretical calculations in Table 4. The 
branching fraction for r~ — > 7r~7r°?7i/^ decay can be predicted by the isospin 
symmetry (CVC) using the experimental results on e"'"e" tt^tt"// [5,22]. 
Our branching fraction agrees with the predictions, B{t^ n^7i^r]Ur) = 
(1.3 ± 0.2) X 10-3 [5] and 1.5 x 10"^ [22]. 

Pich [3], Braaten [23], Li [4] and Aubrecht [24] calculated the branching frac- 
tions of various r decays involving t] meson(s), based on chiral perturbation 
theory, as listed in Table 4. Their results are in fair agreement with our mea- 
surements, particularly taking into account possible uncertainties in theoreti- 
cal predictions. Further detailed studies of the physical dynamics in r decays 
with 7] mesons are required. 

It should be mentioned that the TAUOLA MC program [14] qualitatively 
reproduces our hadronic mass distributions in 7i~7r^riu, t~ K^rju 

and T~ — > K* {892)~ rju decays, as shown in Fig. 7. Further analysis of the 
spectra including various phenomenological models is in progress. 
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Fig. 7. Invariant mass distributions of (a) tt^tt^ and (b) n"!!^?] for r~ — > Ti^TT^rjUT 
decay, (c) K~ri for r~ ^ K^rjUr and (d) Tr^Kgrj for ^ K*~i]i'r, K*~ Kgi:". 
The points with error bars are the data. The normal and fihed histograms indi- 
cate the signal and rr BG MC distributions, respectively. The qq BG is strongly 
suppressed and negligible in our sample. 
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Table 4 

Comparison with theoretical calculations 



Branching fractions 



Mode 



This work 



(1.58 ±0.05 ±0.09) X 10"^ 
(1.35 ±0.03 ±0.07) X 10-3 
(4.6 ±1.1 ±0.4) X 10-^ 
(4.4 ±0.7 ±0.3) X 10"^ 
(1.34 ±0.12 ±0.09) X 10-^ 

< 4.5 X 10^^ 

< 2.5 X 10-^ 

< 3.0 X 10-6 

< 7.4 X 10-6 



Pich [3] Gilman [22] Braaten [23] Li [4] Aubrecht [24] 



r 
r" 

T~ 
T~ 
T~ 



1.2 X 10-^ 

3 X 10-3 1.5 X 10-3 
8.8 X 10-6 
1.1 X 10-^ 

1.6 X lO-'^ 

1.6 X 10-9 
1.1 X 10-9 



2.2 X 10-^ 
1.4 X 10-3 1.9 X 10-3 



1.0 X 10-^ 



1.6 X 10-^ 

7.6 X 10-6 
1.0 X 10-^ 

1.4 X IQ-'^ 

6.9 X 10-9 



